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Abstract 

A new generic dynamical phenomenon of pseudochaos and its relevance 
to the statistical physics both modern as well as traditional one are consid- 
ered and explained in some detail. The pseudochaos is defined as a statistical 
behavior of the dynamical system with discrete energy and/or frequency spec- 
trum. In turn, the statistical behavior is understood as time-reversible but 
nonrecurrent relaxation to some steady state, at average, superimposed with 
irregular fluctuations. The main attention is payed to the most important 
and universal example of pseudochaos, the so-called quantum chaos that is 
dynamical chaos in bounded mesoscopic quantum systems. The quantum 
chaos as a mechanism for implementation of the fundamental correspondence 
principle is also discussed. 

The quantum relaxation localization, a peculiar characteristic implication 
of pseudochaos, is reviewed in both time-dependent and conservative sys- 
tems with special emphasis on the dynamical decoherence of quantum chaotic 
states. Recent results on the peculiar global structure of the energy shell, the 
Green function spectra and the eigenfunctions, both localized and ergodic, in 
a generic conservative quantum system are presented. 

Examples of pseudochaos in classical systems are given including linear 
oscillator and waves, digital computer and completely integrable systems. A 
far-reaching similarity between the dynamics of a few-freedom quantum sys- 
tem at high energy levels (n — > oo) and that of many-freedom one (N — ► oo) 
is also discussed. 



1 A talk given at the Intern. Conference "Nonlinear Dynamics, Chaotic and Complex Systems", 
Zakopane, 1995 
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1 Introduction: the second birth of pseudochaos 



The conception of pseudochaos has been first explicitly introduced [1] in attempts 
to interpret a very controversial phenomenon of quantum chaos, and to understand 
its mechanism and physical meaning. The term itself has been borrowed from the 
theory of the well-known 'pseudorandom number generators' in a digital computer. 
Even though such imitation of the 'true' random quantities is widely used in many 
'numerical experiments', e.g., ones employing the Monte-Carlo techniques, this pseu- 
dorandomness was always considered as a very specific mathematical model of no 
general interest for the fundamental physics. However, in recent numerous attempts 
to understand quantum chaos, which is attracting ever growing attention of many 
researchers (see, e.g., Proceedings of International Conferences [2-4], and a Collec- 
tion of papers [5]), it is becoming more and more clear that this 'specific mechanism' 
provides, in fact, a typical chaotic behavior in physical systems. 

Moreover, from the viewpoint of fundamental physics, the pseudochaos is the 
only kind of chaos principally possible in physical systems of finite dimensions. In 
infinite macroscopic systems of traditional statistical mechanics (TSM) , both classi- 
cal and quantal, particularly, in the principal TSM conception of the thermodinamic 
limit iV — > oo, where N is the number of freedoms, there is no such problem. 
Namely, it has been rigorously proven (see, e.g., Ref. [6]) that, loosely speaking, the 
'true' chaos is a generic phenomenon in this limit even if for any finite N the system 
is completely integrable! 

The discovery of dynamical chaos in finite (and even few-dimensional) classical 
systems - a fundamental breakthrough in recent decades - has crucially changed the 
classical statistical mechanics. By now, this new mechanism for the statistical laws is 
well understood ( but still not very well known), and received the firm mathematical 
foundations in the modern ergodic theory [6]. 

In all success of the latter a 'minor' problem still remains: such a mechanism 
does not work in finite quantum systems that is ones whose motion is bounded in 
the phase space and, hence, whose energy and frequency spectrum is discrete. 

The simplest solution of this problem, which seems to be almost commonly ac- 
cepted currently, is that the dynamical chaos in such systems is simply impossible. 
However, this seemingly obvious 'solution' is, in fact, a trap as it immediately leads 
to a sharp and very profound contradiction with the fundamental correspondence 
principle [7]. We need to choose what to sacrifice, this principle or the 'true' ^clas- 
sical) chaos. I prefer the latter. If the phenomenon of quantum chaos did really 
violate the correspondence principle as some physicists suspect it were, indeed, a 
great discovery since it would mean that the classical mechanics is not the limiting 
case of quantum mechanics but a different separate theory. 'Unfortunately', there 
exists a less radical (but also interesting and important) resolution of this difficulty 
- the pseudochaos - which is the main topic of my talk. 

Within such a philosophical framework the central physical problem is to under- 
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stand the nature and mechanism of dynamical chaos in quantum mechanics. In other 
words, we need the quantum theory of dynamical chaos including the transition to 
the classical limit. Certainly, the quantum chaos is a new dynamical phenomenon 
[7], related but not identical to the classical chaos. We call it pseudochaos, the term 
pseudo intending to emphasize the difference from the 'classical' chaos in the ergodic 
theory. From the physical point of view, I accept here, the latter, strictly speak- 
ing, does not exist in the Nature. So, in the common philosophy of the universal 
quantum mechanics the pseudochaos is the only true dynamical chaos. The classical 
chaos is but a limiting pattern which is, nevertheless, very important both in the 
theory to compare with the real (quantum) chaos and in applications as a very good 
approximation in macroscopic domain as is the whole classical mechanics. Ford calls 
it mathematical chaos as contrasted to the real physical chaos in quantum mechanics 
[10]. 

I emphasize again that the classical chaos is impossible in finite and closed quan- 
tum systems to which my talk is restricted. Particularly, I am not going to discuss 
here the quantum measurement in which, by purpose, the macroscopic (infinite- 
dimensional) processes are involved (see, e.g., Ref.[7]). 

Thus, the physical meaning of the term 'pseudochaos' is principally different 
(and even opposite) to that of 'pseudorandom numbers' in computer. The reason 
for the latter, original, term 'pseudo' was twofold. At the beginning, the first and 
apparently the only meaning was related to the common belief that no dynamical, 
deterministic, system like computer can produce anything random, by definition. 
This delusion has been overcome in the theory of dynamical chaos on the field of 
real numbers. However, the digital computer works on a finite lattice of integers. 
This is qualitatively similar to a quantum behavior [11]. Computer numbers as well 
as quantum variables can be, at most, 'pseudorandom' only as compared to the 
'true' random classical quantities represented by real numbers. But then, a very 
special notion of 'pseudorandom' is sharply scrambling up to the level of a new 
fundamental conception in physics. 

The quantum chaos is a part of quantum dynamics which, in turn, is a particular 
class of dynamical systems. It became a hard physical problem upon discovery and 
understanding of the classical dynamical chaos. To explain the problem I need to 
briefly remind the main peculiarities of the classical chaos, especially those crucial 
in the quantum theory. 

2 Asymptotic chaos in classical mechanics 

There are two equivalent description of classical mechanics or, more generally, of 
any finite-dimensional dynamical systems: via individual trajectories, and via a 
distribution function, or phase-space density for Hamiltonian (most fundamental) 
systems. 
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The trajectory obeys the motion equations, which are generally nonlinear, and 
desribes a particular realization of system's dynamics in dependence on the initial 
conditions. The phase density satisfies the Liouville equation, which is always linear 
whatever the motion equations, and which usually represents the typical (generic) 
dynamical behavior for a given system. Particularly, all zero-measure sets of special 
trajectories are automatically excluded. 

Notice, however, that in some special cases the phase density may display the 
properties absent for trajectories. An interesting example [57] is the correlation 
decay (and, hence, continuous spectrum) for a special initial phase density in a 
completely integrable system. The point is that such a decay is related to the cor- 
relation between different trajectories rather than on the same ones. The trajectory 
spectrum remains descrete, and the corresponding correlation persists. Particularly, 
this explains a surprising phenomenon known as 'echo' which is the revival of such 
correlations upon velocity reversal. It was observed in many cases, and has several 
interesting applications (see, e.g., Ref. [63]). An interesting open question is the ex- 
act conditions for a phase density to represent the trajectory properties which is the 
primary problem in dynamics. 

The strongest statistical properties of a dynamical system are related to the 
local exponential instability of trajectories, as described by the linearized motion 
equations, provided the motion is bounded in phase space. These two conditions 
are sufficient for a rapid mixing of trajectories by the mechanism of 'stretching and 
folding'. For the linear motion equations the combination of both is impossible 
unless the whole phase space of the system is finite. A well-known example of the 
latter is a model described by the linear 'Arnold cat map' [8]: 

p = p + x mod 1 , . 

x = x + p mod 1 

on a unit torus. The motion is exponentially unstable with (positive) Lyapunov's 
exponent A = In [(3 + y/b)/2) > 0, and is bounded by operation mod 1. Notice 
that the linearized motion is described by the same map but without mod 1 that is 
in infinite plane (— oo < dp,dx < oo). It is unbounded and globally unstable but 
perfectly regular, the so-called hyperbolic motion: 

dp = a ■ exp (At) + b ■ exp (—At); dx = c ■ exp (At) + d ■ exp (—At) (2) 

where constants a, b, c, d depend on the initial conditions and on A, and where integer 
t is discrete map's time. Remarkably, the motion (2) is time-reversible but unstable 
in both senses (t — > ±oo). This implies the time reversibility of all the statistical 
properties for the main system (1), a surprising conclusion which is still confusing 
some researchers (see, e.g., Ref. [61]). 

A nontrivial part of the relation between instability and chaos is in that the 
instability must be exponential. A power-law instability is insufficient for chaos. 
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For example, if we change first Eq.(l) to p = p the model becomes completely 
integrable with the oscillation frequency depending on the motion integral p (non- 
linear oscillation). The latter produces linear (in time) instability but the motion 
remains regular (with discrete spectrum). This is a typical property of the com- 
pletely integrable nonlinear oscillations [58] which leads to a confusing difference in 
the dynamical behavior of trajectories and of phase densities as mentioned above. 
Another open question is how to choose the correct time variable for a particular 
dynamical problem [7]. A change of time may convert the exponential instability 
into a power-law one, and vice versa (see, e.g., Ref. [59] for discussion). 

The two above conditions for dynamical chaos can be realized in very simple 
(particularly, few-dimensional) systems like model (1). Another simple example, to 
which I will refer below, is the so-called 'kicked rotator' described by the 'standard 
map' [1,7,9,11]: 

p = p + k ■ sin x; x = x + T ■ p (3) 

also on a torus (x,p mod 2n) or on a cylinder (x mod 2n , — oo < p < oo). This 
model is well studied too, and has many physical applications. The motion on 
cylinder is bounded in one variable only, yet it is sufficient for chaos. 

The exponential instability implies continuous spectrum of the motion which is 
equivalent, loosely speaking, to the mixing, or temporal correlation decay. Appar- 
ently, this is the most important characteristic property in the statistical mechanics 
underlying the principal and universal statistical phenomenon of relaxation to some 
steady state, or statistical equilibrium. 

Aperiodic relaxation is especially clear in the Liouville picture for the phase 
density behavior (see, e.g., Ref. [8]). Consider a certain basis for Liouville's equation, 
for example 

ip mn = exp[2iri(mx + np)\ (4) 

where m,n are any integers in a simple example of model (1). In other words, we 
represent the phase density as a Fourier series: 

f(x, p,t) = J2 Fmn(t)(p mn (x, p) = F mn(0) exp [2n i(m(t)x + n(t)p)\ (5) 

m,n m,n 

Except ipoo any other term of this series has zero total probability, and characterizes 
the spatial correlation in the phase density. Map (1) induces a map for the Fourier 
amplitudes, and for harmonic numbers: 

F = F— n=n+m 
™ mn ' m = m + n [0) 

Remarkably, variables m(t), n(t) obey the same map as for the linearized motion 
equations in variables dx, dp and with the same instability rate A on the infinite 
lattice (m,n). The dynamics of the phase density in the Fourier representation, de- 
scribed by the same Eq.(2) (upon substitution of m, n for dx, dp), is also unbounded, 
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globally unstable, and regular. This is of no surprise as both representations describe 
the local structure of motion. Dynamical chaos is a global phenomenon determined, 
nevertheless, by the microdetails of the initial conditions due to the exponential 
instability of motion [7,13]. Accordingly, in the original phase space the temporal 
density fluctuations are chaotic as are almost all trajectories of map (1). 

The only stationary mode m = n = with the full probability represents in 
this picture the statistical steady state while all the others describe nonstationary 
fluctuations. The latter are another characteristic property of statistical behavior. 
These can be separated from the average statistical relaxation by the so-called 
coarse-graining, or spatial averaging, which is a projection of the phase density on 
a finite (and arbitrarily fine) partition of the phase space. The kinetic (particularly, 
diffusive) description of the statistical relaxation is restricted to such a coarse- 
grained projection only while the fluctuations work as a dynamical generator of 
noise. 

Another elegant method of separating out the average relaxation is the suppres- 
sion of the fluctuations using Prigogine's A operator [61] which provides an invertible 
smoothing of the exact phase density [60] . True, the inverse operator in a nonproper 
one, yet this method could be efficiently used in some theoretical constructions. Con- 
trary to a common belief, it has nothing to do with the time irreversibility [61,62]. 
Moreover, unlike the coarse-grained projection the A-smoothed phase density is as 
reversible as the exact one (principally but not practically, of course). The origin of 
misunderstanding concerning 'irreversibility' is apparently related to the necessary 
restriction on the initial smoothed density which was missed in the theory [62] . Such 
a density is a technical, rather than natural, property of the system, and hence it 
does not need to be arbitrary. A similar operation is often used in quantum me- 
chanics (for different purposes) to convert the Wigner function (the counterpart of 
exact classical phase density) into the so-called Husimi distribution which is the 
expansion in the coherent states (see, e.g., Ref. [7]). 

Nonstationary fluctuations/correlations of the phase density form a stationary 
flow into higher modes |m|, \n\ — > oo (cf. Ref. [12]), and keep the memory of the 
exact initial conditions (see first Eq.(6)) providing time reversibility for the exact 
density. The stationary correlation flow is only possible for the continuous phase 
space which is a characteristic feature of classical mechanics. This allows for the 
asymptotic (t — > ±oo) formulation of the ergodic theory. Notice that both trajec- 
tories and the full density are time-reversible but the latter, unlike the former, is 
nonrecurrent. Reversed relaxation, particularly 'antidiffusion', describes the growth 
of a big fluctuation which is eventually (as t — > — oo) followed by the standard 
relaxation in the opposite direction of time [13]. 
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3 Quantum pseudochaos: a new dimension in the 
ergodic theory 



The dynamical chaos is one limiting case of the modern general theory of dynamical 
systems which describes the statistical properties of deterministic motion (see, e.g., 
Ref. [6]). No doubt, this theory has been developed on the basis of classical me- 
chanics. Yet, as a general mathematical theory, it does not need to be restricted to 
classical mechanics only. Particularly, it can be, and indeed was applied to quantum 
dynamics with a surprising result. Namely, as had been found from the beginning[14] 
and was subsequently well confirmed (see, e.g., Ref.[l,7,ll,22,28]) the quantum me- 
chanics does not typically permit the 'true' (classical-like) chaos. This is because in 
quantum mechanics the energy (and frequency) spectrum of any system, whose mo- 
tion is bounded in phase space, is discrete and its motion is almost periodic. Hence, 
according to the existing ergodic theory, such a quantum dynamics belongs to the 
limiting case of regular motion which is opposite to dynamical chaos. The ultimate 
origin of quantum almost-periodicity is in the discreteness of the phase space itself 
(or, in a more formal language, in the noncommutative geometry of the latter) which 
is at the basis of quantum physics and directly related to the fundamental uncer- 
tainty principle. Yet, another fundamental principle, the correspondence principle, 
requires the transition to classical mechanics in all cases including the dynamical 
chaos with all its peculiar properties. 

Now, the principal question to be answered reads: where is the expected quantum 
chaos in the ergodic theory? The answer to this question [7,11,13] (not commonly 
accepted as yet) was concluded from a simple observation (principally well-known 
but never comprehended enough) that the sharp border between the discrete and 
continuous spectrum is physically meaningful in the limit \t\ — > oo only, the condition 
actually assumed in the ergodic theory. Hence, to understand the quantum chaos the 
existing ergodic theory needs some modification by introducing a new 'dimension', 
the time. In other words, a new and central problem in the ergodic theory becomes 
the finite-time statistical properties of a dynamical system, both quantal as well as 
classical. 

Within a finite time the discrete spectrum is dynamically equivalent to the con- 
tinuous one, thus providing much stronger statistical properties of the motion than 
it was (and still is) expected in the ergodic theory in case of discrete spectrum. 
It turns out that the motion with discrete spectrum may exhibit all the statistical 
properties of the classical chaos but only on some finite time scales. 

The absence of the classical-like chaos in quantum dynamics apparently con- 
tradicts not only the correspondence principle but also the fundamental statistical 
nature of quantum mechanics. However, even though the random element in quan- 
tum mechanics ('quantum jumps') is unavoidable, indeed, it can be singled out and 
separated from the proper quantum processes. Namely, the fundamental random- 
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ness in quantum mechanics is related only to a very specific event - the quantum 
measurement - which, in a sense, is foreign to the proper quantum system itself. 
This allows to divide the whole problem of quantum dynamics in two qualitatively 
different parts: 

• The proper quantum dynamics as described by a very specific dynamical vari- 
able, the wavefunction ip(t) obeying some deterministic equation, for example, 
the Schrodinger equation. The discussion below will be limited to this part 
only. 

• The quantum measurement including the registration of the result and, hence, 
the collapse of the ip function which still remains a very vague issue to the 
extent that there is no common agreement even on the question whether this 
is a real physical problem or an ill-posed one, so that the Copenhagen inter- 
pretation of quantum mechanics answers all the 'admissible' questions. In any 
event, there exists as yet no dynamical description of the quantum measure- 
ment including the ip collapse. 

Recent breakthrough in the understanding of quantum chaos has been achieved, 
particularly, due to the above philosophy of separating out the dynamical part of 
quantum mechanics. Such a philosophy is accepted, explicitly or more often implic- 
itly, by most researchers in this field. 

3.1 Time scales of pseudochaos 

The existing ergodic theory is asymptotic in time, and thus has no explicit time scales 
at all 0. There are two reasons for this. One is technical: it is much simpler to derive 
the asymptotic relations than to obtain rigorous finite-time estimates. Another rea- 
son is more profound. All statements in the ergodic theory hold true up to measure 
zero that is excluding some peculiar nongeneric sets of zero measure. Even this min- 
imal imperfection of the theory did not seem completely satisfactory but has been 
'swallowed' eventually and is now commonly tolerated even among mathematicians, 
let alone physicists. In a finite-time theory all these exceptions acquire a small but 
finite measure which would be apparently 'unbearable' (for mathematicians). Yet, 
there is a standard mathematical 'trick' for avoiding both these difficulties. 

The most important time scale t R in quantum chaos is given by the general 
estimate [7,11]: 

Q a 

hi(ut R ) ~ InQ, t R ~ — ~ p < Ph (7) 

2 Asymptotic statements in the ergodic theory should not be always understood literally to avoid 
physical misconceptions (see, e.g., Addendum in second Ref.[12]). Actually, the classical chaos has 
also time scales, for example, a dynamical one (~ A -1 ) [7]. 
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where u and a ~ 1 are system-dependent parameters, and Q > 1 stands for some 
big (in semiclassical region) quantum parameter. It may be, e.g., a quantum number 
Q — I/H related to a characteristic action variable I or the total number of states for 
the bounded quantum motion in a phase space domain of volume r: Q ~ T/(2tt) n . 
Here and below I set fi — 1. 

This scale is called the relaxation time scale referring to one of the principal 
properties of the chaos - statistical relaxation to some steady state. The physical 
meaning of this scale is principally simple, and is directly related to the fundamental 
uncertainty principle (At- AE ~ 1) as implemented in the second Eq.(7) where pn 
is the full average energy level density (also called Heisenberg time). For t£tn the 
discrete spectrum is not resolved, and the statistical relaxation follows the classical 
(limiting) behavior. This is just the 'gap' in the ergodic theory (supplemented with 
the additional dimension, the time) where the pseudochaos, particularly quantum 
chaos, dwells. A more accurate estimate relates i# to a part po of the level density. 
This is the density of the so-called operative eigenstates only, that is those which 
are actually present in a particular quantum state ip, and which do actually control 
its dynamics. 

The formal trick mentioned above is to consider not finite-time relations, we 
really need in physics, but rather the special conditional limit: 

t, Q -> oo , t r = = const (8) 

where tr is a new dimensionless time. The double limit (8) (unlike the single one 
Q — > oo) is not the classical mechanics which holds true, in this representation, for 
tr ^ 1 and with respect to the statistical relaxation only. For tr £ 1 the behavior 
becomes essentially quantum (even in the limit Q — > oo !) and is called nowadays 
mesoscopics. Particularly, the quantum steady state is generally quite different from 
the classical statistical equilibrium in that the former may be localized (under certain 
conditions) that is nonergodic in spite of classical ergodicity. 

Another important difference is in fluctuations which are also a characteristic 
property of chaotic behavior. In comparison with classical mechanics the quantum 
ip(t) plays, in this respect, an intermediate role between the classical trajectory with 
big relative fluctuations ~ 1 and the coarse-grained classical phase density with no 
fluctuations at all. Unlike both the fluctuations of ip(t), or rather those of averages 
in a quantum state ip(t), are typically ~ d H where du is the number of operative 
eigenstates associated with quantum state ip which the former may be also called the 
Hilbert dimension of state ip. In other words, chaotic ip(t) represents statistically a 
finite ensemble of ~ du independent systems even though formally ip(t) describes a 
single system. The fluctuations clearly demonstrate the difference between physical 
time t and auxiliary variable r: in the double limit (t, Q — > oo) the fluctuations 
vanish, and one needs a new 'trick' to recover them for a finite Q. 

The relaxation time scale should not be confused with the Poincare recurrence 
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time tp ^ tp which is typically much longer, and which sharply increases with 
decreasing the recurrence domain. Time scale tp characterizes big fluctuations (for 
both the classical trajectory, but not the phase density, and for the quantum if)) 
of which recurrences is a particular case. Unlike this, tp characterizes the average 
relaxation process. Rare recurrences, the more rare the larger quantum parameter 
Q, make quantum relaxation similar to the classical nonrecurrent one. 

More strong statistical properties than relaxation and fluctuations are related 
in the ergodic theory to the exponential instability of motion. The importance 
of those stronger properties for the statistical mechanics is not completely clear 
[56]. Nevertheless, in accordance with the correspondence principle, those stronger 
properties are also present in quantum chaos as well but on a much shorter time 

At r ~ InQ (9) 

where A is the classical Lyapunov exponent. This time scale was discovered and 
partly explained in Ref.[15] (see also Ref. [7,11]). We call it random time scale. 
Indeed, according to the Ehrenfest theorem the motion of a narrow wave packet 
follows the beam of classical trajectories as long as the packet remains narrow, and 
hence it is as random as in the classical limit. Even though the random time scale 
is very short, it grows indefinitely as Q — > oo. Thus, a temporary, finite-time 
quantum pseudochaos turns into the classical dynamical chaos in accordance with 
the correspondence principle. 

Again, we may consider the conditional limit: 

t,Q—>oo, r r = j = const (10) 

Notice that new scaled time r r is different from the previous one Tp in Eq.(8). 

Particularly, if we fix time t, then in the limit Q — > oo we obtain the transition 
to the classical instability in accordance with the correspondence principle while for 
Q fixed, and t — > oo we have the proper quantum evolution in time. For example, 
the quantum Lyapunov exponent 

The quantum instability (A 9 > 0) was observed in numerical experiments [7,16]. 
What does terminate the instability for t £ t r ? A simple explanation is suggested by 
the classical picture of the phase density evolution on the integer Fourier lattice m, n 
discussed above for model (1). Classical Fourier harmonics m, n are of a kinematical 
nature without any a priori dynamical restriction. Particularly, they can go, and do 
so for a chaotic motion, arbitrarily large which corresponds to a continuous classical 
phase space. On the contrary, the quantum phase space is discrete. At first glance, 
the quantum wave packet stretching/squeezing, similar to the classical one, does 
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not seem to be principally restricted since only 2-dimensional area (per freedom) is 
bounded in quantum mechanics. However, Fourier harmonics of the quantum phase 
density (Wigner function) are directly related to the quantum dynamical variables, 
particularly, to the action variables whose values are restricted by the quantum 
parameter Q, hence estimate (9). In a simple model (1) this is related to a finite size 
of the whole phase space. Generally, in a conservative system with even infinite phase 
space the restriction is imposed by the energy conservation. Numerical experiments 
reveal that the original wave packet, after a considerable stretching similar to the 
classical one, is rapidly destroyed. Namely, it gets split into many new small packets 
[7,16]. The mechanism of this sharp 'disrupture' of the classical-like motion is not 
completely clear (for a possible explanation see Ref.[7,17]). The resulting picture is 
qualitatively similar to that for the classical phase density, the main difference being 
in the spatial fluctuation scale bounded now from below by 1/Q. Nevertheless, the 
quantum phase density can be also decomposed into a coarse-grained average part, 
and the fluctuations. An important implication of this picture for the wave packet 
time evolution is the rapid and complete destruction of the so-called generalized 
coherent states [18] in quantum chaos. 

In quasiclassical region (Q > 1) the scale t r <C t R . This leads to a surprising 
conclusion that the quantum diffusion and relaxation are dynamically stable con- 
trary to the classical behavior. It suggests, in turn, that the motion instability is, 
generally, not important during statistical relaxation. However, the foregoing cor- 
relation decay on the short time scale t r is crucial for the statistical properties of 
quantum dynamics. 

Dynamical stability of quantum diffusion has been proved in striking numerical 
experiments with time reversal [19]. In a classical chaotic system the diffusion is 
immediately recovered due to numerical 'errors' (not random !) amplified by the 
local instability. On the contrary, the quantum 'antidiffusion' proceeds until the 
system passes, to a very high accuracy, the initial state, and only then the normal 
diffusion is restored. The stability of quantun chaos on relaxation time scale is 
comprehensible as the random time scale is much shorter. Yet, the accuracy of 
the reversal (up to ~ 1CT 15 (!) ) is surprising. Apparently, this is explained by a 
relatively large size of the quantum wave packet as compared to the unavoidable 
rounding-off errors unlike the classical computer trajectory which is just of that size 
[20]. In the standard map (3) (upon quantization) the size of the optimal, least- 
spreading, wave packet Ax ~ VT [11]. On the other hand, any quantity in the 
computer must well exceed the rounding-off error 5< 1. Particularly, T 5, and 
(Axf/5 2 >(T/5)5- 1 ^>l. 

3.2 Classical— like relaxation and residual fluctuations 

The relaxation time scale tn is the most important of the two considered above for 
two reasons. First, it is much longer than t r , and second, it is related to the principal 
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process of statistical relaxation which is the basis of statistical mechanics. The short 
scale t r was interpreted in Ref.[15] (see also Ref.[26]) as a limit for the classical-like 
behavior of chaotic quantum motion. Subsequently, it was found that the method of 
quasiclassical quantization can be extended on a much longer time [11,27]. However, 
the physics on both time scales is qualitatively different: dynamical instability on 
scale t r , and statistical relaxation afterwards. 

On the whole scale tu the discrete pseudochaos spectrum is not resolved, and 
relaxation follows the classical law. Consider, for example, model (3), the standard 
map on a torus with total number of quantum states Q, and p, x as the action-angle 
variables. 

If perturbation parameter k^Q the relaxation to ergodic steady state in this 
model, as well as in model (1), is very quick, with characteristic relaxation time 
t e ~ 1 (iterations). Such regime does often take place in many physical systems. 
Here I consider another, more interesting for the problem of pseudochaos, case, 
namely the diffusive relaxation which occurs for a sufficiently weak perturbation 

k < Q (12) 

In the classical limit this relaxation is descibed by the standard diffusion equation 

where f(p,t) =< f(p,x,t) > x is a coarse-grained phase density (averaged over x), 
and 

D = < (A " )2 > » k - (14) 
~b 2 

the diffusion rate. The latter expression holds for the standard map if K = kT ^> 1 
which is also the condition for the global chaos in this model [9]. The relaxation to 
the ergodic steady state f s — 1/Q is exponential with characteristic time 

t - Q2 ~ Q2 (15) 

In diffusive regime (k <C Q) this time t e 3> 1. That average relaxation is stable and 
regular in spite of underlying chaotic dynamics. 

The quantized standard map ip = Uip is described by a unitary operator 

( Tp 2 \ 

U = exp I — i — ^— j • exp { — ik- cosx) (16) 

on a cylinder (Q — > oo) [14], where p = —id/dx, and by a similar but somewhat 
more complicated expression on a torus [21]. 

There are three quantum parameters in this model: perturbation k, period T 
and size Q but only two classical combinations remain: perturbation K — k-T, and 
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classical size M = TQ/2-k which is the number of classical resonances over the torus. 
Notice that the quantum dynamics is generally more rich than the classical one as 
the former depends on an extra parameter. It is, of course, another representation 
of Planck's constant which I have set H — 1. This is why in the quantized standard 
map we need both parameters, k and T, separately and cannot combine them in a 
single classical parameter K. 

The quasiclassical region, where we expect quantum chaos, corresponds to T — > 
0, k — > oo, Q — > oo while the classical parameters K = const and M = const are 
fixed. 

A technical difficulty in evaluating t# for a particular dynamical problem is in 
that the density p depends, in turn, on the dynamics. So, we have to solve a 
self-consistent problem. For the standard map the answer is known (see Ref.[7]): 



This is a remarkable relation as it connects essentially quantum characteristics 
(tR, Po) with the classical diffusion rate D (14). 

The quantum diffusion rate depends on the scaled (dimensionless) time tr (8), 
and is given by 



This is an example of scaling in discrete spectrum which eventually stops the quan- 
tum diffusion. 

The character of the steady state crucially depends on the ratio tnjt e . Define 
the ergodicity parameter A as [7] 



Consider, first, the case A> 1 when the time scale tn is long enough to allow for 
the completion of the classical-like relaxation. In this case the final steady state as 
well as all the eigenfunctions are ergodic that is the corresponding Wigner functions 
are close to the classical microcanonical distribution in phase space. This region is 
inevitably reached if the classical parameter K/M is kept fixed while the quantum 
parameter k — > oo in agreement with the Shnirelman theorem or, better to say, with 
a physical generalization of that [23] . It is called the far quasiclassical asymptotics. 

The principal difference of the quantum ergodic state from the classical one is 
residual fluctuations in the former. In quasiclassical region the chaotic quantum 
steady state is a superposition of very many eigenfunctions. As a result almost any 
physical quantity fluctuates in time. Even in discrete spectrum we are considering 
here these fluctuations are very irregular. In case of a classical-like ergodic steady 
state all Q eigenfunctions essentially contribute to the fluctuations. Moreover, we 



Po 



2D 



(17) 




(18) 




(19) 
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would expect their contributions to be statistically almost independent. Hence, the 
fluctuations should scale ~ Q^ 1 ^ 2 = d^ 2 where du is the Hilbert dimension of the 
ergodic state. This is the case, indeed, according to numerical experiments [29]. For 
example, the energy fluctuations were found to follow a simple relation 



E s ' v/Q 

where 



(20) 



< p 2 > 



(AE S ) 2 = &(t) - E 2 ; E s = E(t)- E(t) = (21) 

Here the bar indicates time averaging over a sufficiently long time interval (>> t e ), 
and the brackets denote usual average over the quantum state. 

Dependence (20) suggests the complete quantum decoherence in the final steady 
state for any initial state even though the steady state is formally a pure quantum 
state. For Q ^> 1 the fluctuations are small, so that statistically the quantum 
relaxation is nonrecurrent. The decoherence of a chaotic quantum state is also 
confirmed by the independence (up to small fluctuations) of the final steady-state 
energy E s on the initial E(0). Since any particular initial quantum state is strongly 
coherent the decoherence is a result of the quantum chaos. It is called dynamical 
decoherence. This is one of the most important results of the studies in quantum 
chaos. 



3.3 Mesoscopics: quantum behavior in quasiclassical region 

If ergodicity parameter A <C 1 is small all the eigenstates and the steady state are 
non-ergodic, or localized. This is because the scale tn is not long enough to support 
the classical-like diffusion which stops before the classical relaxation is completed. 
For this reason it is also called the quantum diffusion localization. As a result the 
structure of eigenfunctions and of the steady state remains essentially quantum, no 
matter how large is the quantum parameter k — > oo. This is called intermediate 
quasiclassical asymptotics or mesoscopic domain. Particularly, it corresponds to 
K > 1 fixed, k — > oo and M — > oo while A 1 remains small. 

The popular term 'mesoscopic' means here some intermediate behavior between 
classical and quantum one. In other words, in mesoscopic phenomenona both classi- 
cal and quantum features are combined simulteneously. Again, the correspondence 
principle requires transition to the completely classical behavior. This is, indeed, the 
case as the mesoscopic phenomena occur in the region where quantum parameter 
k ^> 1 is already very big but still less than a certain critical value (corresponding 
to A ~ 1) which determines the border of transition to the fully classical behavior 
(far quasiclassical asymptotics). 

If A 1 is very small the shape of the localized eigenstates is asymptotically 
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exponential [7], and can be approximately described by a simple expression [41]: 

/.,.M^<i^( P )i 2 >^ coBh|2( 2 /:' pm)/il (22) 

The localized steady state has a similar but somewhat more complicated shape 
[11,31]. This is a simple approximation superimposed with big fluctuations. The 
parameter / is called localization length. Interestingly, the two localization lengths 
(l s for steady state and I for eigenfunctions) are rather different [11]: 

l s ps D while / ps — (23) 

2 

because of big fluctuations. 

In terms of localization length the region of mesoscopic phenomena is defined by 
the double inequality: 

1 « I < Q (24) 

The left inequality is a classical feature of the state while the right one refers to 
quantum effects. The combination of both allows, particularly, for a classical de- 
scription, at least in the standard map, of statistical relaxation to the quantum 
steady state by a phenomenological diffusion equation [7,24] for the Green function: 

dg(v, a) 1 d 2 g . dg /n . 

Here g(u, 0) = 0)| 2 = 5(u — is ) and 

V= 2D" cr = ln ( 1 + ^)> TR = ^3 ^ 
The additional drift term in the diffusion equation with 

B(y) ps sign(z/ - u ) = ±1 (27) 

describes the so-called quantum coherent backscattering which is the dynamical 
mechanism of localization. 

The solution of Eq.(25) reads [7]: 



9iy, a ) = —;= ex P 



(A + a) 



+ exp ( - 4A) • erfc ( ^=^ J ( 28 ) 



where A — \u — uq\. 

Asymptotically, as a — > oo, the Green function g(u, a) — > 2 exp (—4 A) = g s 
approaches the localized steady state g s , exponentially in o but only as a power-law 
in physical time r R or t (g—g s ~ 1/ t r ) . This is the effect of discrete motion spectrum. 
Numerical experiments confirm prediction (28), at least, to the logarithmic accuracy 
~ a ps hi 77? [7,22]. 
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A physical example of localization is the quantum suppression of diffusive photo- 
effect in Hydrogen atom [25]. Depending on parameters the suppression may occur 
no matter how large the atomic quantum numbers. This is a typical mesoscopic 
phenomenon which had been predicted by the theory of quantum chaos, and was 
subsequently observed in laboratory experiments. 

One might expect that in case of localization (D <C Q) the fluctuations would 
scale like Z -1 / 2 ~ k~ l as the number of eigenfunctions coupled in the localized steady 
state is ~ /. This is, however, not the case as was found already in first numerical 
experiments (see Ref.[7]). According to more accurate data [29] the fluctuations are 
described by the relation 

^ = A - (29 ) 

with fitting parameters 7 = 0.55, a = 0.65. For a nonergodic state the Hilbert 
dimension can be defined as (see Ref.[7]) 

*H=\-j \^iP)\' d P = J f(P)dP (3°) 

where f(p) is smoothed (coarse-grained) density, and the factor 1/3 accouts for 
ip fluctuations [30]. In case of exponential localization (22) dn ~ 7r 2 //4. The 
most important parameter 7 here is about twice as small compared to the expected 
value 7 = 1. This result suggests some fractal properties of localized eigenfunctions 
and/or of their spectra. To put it another way, a slow fluctuation decay (29) implies 
incomplete quantum decoherence which can be characterized by the number d s of 
statistically independent components in the steady state [7]. Then, from Eqs.(20) 
and (29) we obtain in the two limits: 

T H ~\ d Jy = 2.4 A « 1 (31) 

This result was confirmed in Ref.[31] for a band random matrix model. 

The phenomenon of quantum diffusion localization explains also the limitation 
of quantum instability in systems with infinite phase space like the standard map 
on a cylinder. Indeed, the maximal number of coupled states here is determined by 
the localization length whatever the total number of states in the system. Hence, we 
should substitute quantum parameter Q ~ / ~ k 2 in estimate (9). Even if localiza- 
tion does not take place (e.g., for standard map with parameter k(p) depending on 
p, see Ref. [7,11]), so that the quantum diffusion doesn't stop at all and the quantum 
spectrum becomes continuous, the number of coupled states increases with time as 
a power law only (Ap ~ y/i), and hence the quantum Lyapunov exponent A q — > 
vanishes on the relaxation time scale. Only if the action variables grow exponentially 
the instability rate A q remains finite, and the quantum chaos becomes asymptotical 
like in the classical limit (see Ref. [11,32] for such 'exotic' models). 
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3.4 Examples of pseudochaos in classical mechanics 

The pseudochaos is a new generic dynamical phenomenon missed in the ergodic 
theory. No doubt, the most important particular case of pseudochaos is the quantum 
chaos. Nevertheless, pseudochaos occurs in classical mechanics as well. Here are a 
few examples of classical pseudochaos which may help to understand the physical 
nature of quantum chaos. Besides, it unveils new features of classical dynamics as 
well. 

Linear waves is the most close to quantum mechanics example of pseudochaos 
(see, e.g., Ref. [33]). I remind that here only a part of quantum dynamics is discussed, 
one described, e.g., by the Schrodinger equation which is a linear wave equation. 
For this reason the quantum chaos is called sometime wave chaos [34]. Classical 
electromagnetic waves are used in laboratory experiments as a physical model for 
quantum chaos [35]. The 'classical' limit corresponds here to the geometrical optics, 
and the 'quantum' parameter Q = L/X is the ratio of a characteristic size L of the 
system to wave length A. As is well known in optics, no matter how large is the ratio 
L/X the diffraction patterns prevail at a sufficiently far distance R^L 2 /X. This is 
a sort of relaxation scale: R/X ~ Q 2 . 

Linear oscillator (many-dimensional) is also a particular representation of 
waves. A broad class of quantum systems can be reduced to this model [36] . Statisti- 
cal properties of linear oscillator, particularly in the thermodynamic limit (N — > oo), 
were studied in Ref. [37] in the framework of TSM. On the other hand, the theory of 
quantum chaos suggests more rich behavior for a big but finite N, particularly, the 
characteristic time scales for the harmonic oscillations [38], the number of freedoms 
N playing a role of the 'quantum' parameter. 

Completely integrable nonlinear systems also reveal pseudochaotic behav- 
ior. An example of statistical relaxation in the Toda lattice had been presented in 
Ref. [39] much before the problem of quantum chaos arose. Moreover, the strongest 
statistical properties in the limit N — > oo, including one equivalent to the expo- 
nential instability (the so-called if-property) were rigorously proved just for the 
systems completely integrable for any finite N (see Ref. [6]). 

Digital computer is a very specific classical dynamical system whose prop- 
erties are extremely important in view of the ever growing interest to numerical 
experiments covering now all branches of science and beyond. The computer is the 
'overquantized' system in that any quantity here is discrete while in quantum me- 
chanics only the product of two conjugated variables does so. 'Quantum' parameter 
here Q = M which is the largest computer integer, and the short time scale (9) 
t r ~ InM, the number of digits in the computer word [11]. Owing to the discrete- 
ness, any dynamical trajectory in computer becomes eventually periodic, the effect 
well known in the theory and practice of pseudorandom number generators. One 
should take all necessary precautions to exclude such computer artifacts in numer- 
ical experiments (see, e.g., Ref.[40,64]). On the mathematical part, the periodic 



17 



approximations in dynamical systems are also considered in the ergodic theory, ap- 
parently without any relation to pseudochaos in quantum mechanics or computer 
[6]. 

The computer pseudochaos seems to me most convincing argument for the re- 
searchers who are still reluctant to accept the quantum chaos as, at least, a kind of 
chaos, insisting that only the classical-like (asymptotic) chaos deserves this name, 
the same chaos which was (and is) studied to a large extent just on computer, that 
is the chaos inferred from a pseudochaos! 

4 Statistical theory of pseudochaos: 
random matrices 

The complete solution of the dynamical quantum problem can be obtained via di- 
agonalization of the Hamiltonian to find the energy (or quasi-energy) eigenvalues 
and eigenfunctions. The evolution of any quantity is, then, expressed as a sum over 
the eigenfunctions. For example, the energy time dependence is 

E(t) = °rn c* m , E mm , exp [i (E rn - E m ,) t] (32) 

mm' 

where E mm i are the matrix elements, and the initial state in momentum representa- 
tion is ip(n, 0) = c m <fm(n)- For chaotic motion the dependence is generally very 

m 

complicated but the statistical properties of the motion can be inferred from the 
statistics of eigenfunctions <p m (n) (and hence of the matrix elements E mm r), and of 
eigenvalues E m . 

By now, there exists a well-developed random matrix theory (RMT, see, e.g., 
Ref. [43]) which describes the average properties of a typical quantum system with 
a given symmetry of the Hamiltonian. At the beginning the object of this theory 
was assumed to be a very complicated, particularly many-dimensional, quantum 
system as a representative of a certain statistical ensemble. With understanding 
the phenomenon of dynamical chaos it became clear that the number of system's 
freedoms is irrelevant. Instead, the number of quantum states (quantum parameter 
Q) is of importance provided the dynamical chaos in the classical limit. 

This approach to the theory of complex quantum systems like atomic nuclei had 
been taken by Wigner [42] 40 years ago, much before the problem of quantum chaos 
was realized. He introduced the so-called band random matrices (BRM) which were 
most suitable to account for the structure of conservative systems. However, due to 
severe mathematical difficulties, RMT immediately turned to a much simpler case 
of statistically homogeneous (full) matrices, for which impressive theoretical results 
have been achieved [43]. The price was that the full matrices describe the local 
chaotic structure only, the limitation especially inacceptable in atoms [30,44]. Only 
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recently the interest of some researchers turned back to the original Wigner BRM 
[45,48]. 

One of the main results in the studies of quantum chaos was the discovery of 
quantum diffusion localization as a mesoscopic quasiclassical phenomenon. This phe- 
nomenon, discussed above, has been well studied and confirmed by many researchers 
for the dynamical models described by maps. Contrary to a common belief the maps 
describe not only time-dependent systems but also conservative ones in the form of 
Poincare maps (see, e.g., Ref. [49]). Nevertheless, to my knowledge no direct studies 
of the quantum localization in conservative systems have been undertaken as yet, 
either in laboratory or even in numerical experiments. Moreover, the very existence 
of quantum localization in conservative systems is challenged [50] . Here, I briefly de- 
scribe recent results concerning the structure of the localized quantum chaos in the 
momentum space of a generic few-freedom conservative system which is classically 
strongly chaotic, particularly, ergodic on a compact energy surface [41]. 

Generally, RMT is a statistical theory of the systems with discrete energy spec- 
trum. This is just the principal property of the quantum dynamical chaos [7]. Thus, 
RMT turned out to become accidentally (!) a statistical theory for the coming quan- 
tum chaos. Remarkably, this statistical theory does not include any time-dependent 
noise that is any coupling to a thermal bath, the standard element of most statisti- 
cal theories. Moreover, a single matrix from a given statistical ensemble represents 
the typical (generic) dynamical system of a given class characterized by the matrix 
parameters. This makes an important bridge between dynamical and statistical de- 
scription of the quantum chaos. In matrix representation the similarity between the 
problem of quantum diffusion localization in momentum space and the well-known 
dual problem of Anderson localization in configurational space of disordered solids 
[53,54] is especially clear and instructive. 

Consider real Hamiltonian matrices of a rather general type 

H mn = H nn 5 mn + v mn m,n = 1, ...,JV (33) 

where off-diagonal matrix elements v mn = v nm are random and statistically inde- 
pendent with < v mn >= and < v 2 mn >= v 2 for \m — n\ < b, and are zero otherwise. 
The most important characteristic of these Wigner band random matrices (WBRM) 
is the average energy level density p defined by the relation 

- = {Hmm — H m i m i) (34) 

p 

where m' — to — 1. The averaging here and below is understood either over disorder 
that is over many random matrices or within a single sufficiently large matrix. Both 
are equivalent owing to the assumed statistical independence of matrix elements. 
In other words, many matrices are statistically equivalent to a big one. Quantum 
numbers to, n are generally arbitrary but we will have in mind ones related to the 
action variables, thus considering the quantum structure in the momentum space. 
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The basis in which the matrix elements are calculated is usually assumed to cor- 
respond to a completely integrable system with iV quantum numbers where N is 
the number of freedoms. By ordering the basis states in energy we can represent 
N quantum numbers by a single one related to the energy which is also an action 
variable. 

In the classical limit the definition of WBRM (33) corresponds to the standard 
Hamiltonian H = Hq + V where the perturbation V is usually assumed to be 
sufficiently small while the unperturbed Hamiltonian Hq is completely integrable. 

Quantum model (33) is defined by 3 independent physical parameters above: 
p, v, and b. The fourth parameter, matrix size Q, is considered to be technical in 
this model provided Q 3> d e (see Eq. (38) below) is big enough to avoid the boundary 
effects. 

In terms of unperturbed energy E the classical chaotic trajectory of a given 
total energy E = const fills up the energy shell AE = AV with the ergodic (micro- 
canonical) measure w e depending on a particular perturbation function V. In the 
quantum system this measure characterizes the shape (distribution) of the ergodic 
eigenfunction (EF) in the unperturbed basis. Conversly, if we keep fixed the unper- 
turbed energy E = const, the measure w e describes the band of energy surfaces 
E = const whose trajectories reach the unperturbed energy Eq. In a quantum sys- 
tem the measure w e in the latter case corresponds to the energy spectrum of a Green 
function (GS) with initial energy E . This characteristic was originally introduced 
also by Wigner [42] as the 'strength function', the term still in use in nuclear physics. 
Now it is called also the 'local density of (eigen) states'. 

For a typical perturbation, represented by WBRM, w e depends on the Wigner 
parameter [42] q = (pv) 2 /b. In the two limits [42] (see also Ref.[7,47]) 



W e (E) 



4^ V^lc - E\ \E\ < E sc , q » 1 



' 2-arctan(l/7rg)' 1^1 - E BW, Q < 1 



(35) 



provided n = pv £ 1 which is the condition for coupling neighboring unperturbed 
states by the perturbation. In opposite case n <^ 1 the impact of perturbation 
is negligible which is called perturbative localization. The latter is a well known 
quantum effect but not one we are interested in (for chaotic phenomena it was first 
considered in Ref.[51]). What is less known that for the coupling of all unperturbed 
states within the Hamiltonian band a stronger condition is required, namely: 

?7 £ Vb, or q £ 1 (36) 

This is a simple estimate in the first order of perturbation theory. Indeed, the 
coupling is ~ V/5E. Within the band the energy detuning 5E ~ b/p while the total 
random perturbation V ~ vy/b, hence estimate (36). In opposite case q £ 1 a partial 
perturbative localization takes place which is also a quantum phenomenon, and again 
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not one we have in mind speaking about the quantum localization. The mechanism 
of perturbative localization is relatively simple and straightforward. This quantum 
effect is completely absent only in the first, semicircle (SC), limit of Eq.(35) where 
the width of the energy shell AE = 2E SC = 2VSbv 2 = A^/2qE b > E b , and E b = b/p 
is the half-width of the Hamiltonian matrix band in energy. The latter inequality 
allows for diffusive quantum motion within the energy shell as a single random jump 
is ~ b <C pAE. The quantum localization under consideration here is related just 
to the localization (suppression) of quantum diffusion by the interference effects in 
discrete spectrum (see, e.g., Ref. [7]). Notice that the SC width immediately follows 
from the above estimate for perturbative localization: SE ~ AE ~ v\/b. 

In the second, Breit - Wigner (BW), limit of Eq.(35) the full size of the energy 
shell AE = 2E B w = ZE b is equal to that of the Hamiltonian band. However, due to 
the partial perturbative localization explained above the main peak of the quantum 
ergodic measure is considerably more narrow, with the width T = 2-npv 2 = 2-xqE b <C 
E b . This is again in accordance with the same simple estimate: 5E ~ T ~ v^fpT. 

To the best of my knowledge, the quantum distributions (35) were theoretically 
derived and studied for GSa only. Classically, the measure w e seems to be the same 
for both E = const and E = const as determined by the same perturbation V. One 
of the main recent results [41] in the studies of WBRM is that the classical symmetry 
between EFs (E = const) and GSa (E = const) is generally lost in quantum 
mechanics. Namely, in ergodic case such a statistical symmetry still persists, yet the 
quantum localization drastically violates the symmetry producing a very intricate 
and unusual global structure of the quantum chaos. 

In a sense, the conservative system is always localized (finite AE) even for ergodic 
motion. This is the origin of misunderstanding sometimes (see, e.g., Ref. [52]). In 
fact, such a classical localization is a trivial consequence of energy conservation as 
was explained above. It persists, of course, in the classical limit as well. Here we 
are interested in the quantum localization explained above called simply localization 
below. 

Similar to maps the localization in conservative systems also depends on the 
ergodicity parameter (cf. Eq.(19)): 

b 2 ab 3 / 2 
A = a — = — -= — 

d e A\j2cn 

Here the ergodicity corresponds not to the total number of states Q as in maps (19) 
but to that within the energy shell of width AE: 

d e = c- p{AE) sc = 4:cnV2b (38) 

Hilbert dimension d e is also called the ergodic localization length as a measure of the 
maximal number of basis states (BS) coupled by the perturbation in case of ergodic 
motion. Numerical factor c ~ 0.92 is directly calculated from the limiting expression 
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(37) 



(35) for a particular definition of d in Eq.(30). Relation (38) is valid formally in 
the SC region (q ^> 1) only but, according to computations, has still the accuracy 
within a few per cent down to q ~ 0.4. 

Parameter A (37) had been found in Ref.[48] and implicitely used there (without 
any relation to ergodicity). It was explained in details in Ref.[45] where factor 
a w 1.2 was also calculated numerically. 

Localization is characterized by the parameter 

(3 d = » 1 - e- A < 1 (39) 

Here d stands for the actual average localization length of EFs measured according 
to the same definition (30). Empirical relation (39) has been found in numerical 
experiments [45] to hold in the whole interval A < 2.5 there, and was confirmed in 
Ref.[41] up to A « 7. 

In the BW region d e = npT = 2n 2 bq, and A ^ ab/(2n 2 q) > 1 as q < 1 (35) 
and b ^> 1 in quasiclassics. Hence, localization is only possible in the SC domain 
which was studied in Ref.[41]. 

The numerical results [41] were obtained from two individual matrices: the main 
one for the localized case with parameters 

A = 0.23; q = 90; Q = 2400; v = 0.1; b = 10; p = 300; r] = 30; d e = 500 

and additional one for ergodic case with parameters 

A = 3.6; q = 1; Q = 2560; v = 0.1; b = 16; p = 40; r] = 4; d e = 84 

All results are entirely contained in the EF matrix c mn which relates the eigen- 
functions ip m to unperturbed basis states ip n : 

i>m = C mn-fn, W m {n) = |^(n)| 2 = C^ n = W mn (40) 

n 

in momentum representation, and in the corresponding string of eigenvalues E m m 
m/p. From the matrix c mn the statistics of both EFs as well as GSa was evaluated. In 
order to suppress big fluctuations in individual distributions the averaging over 300 
of them in the central part of the matrix was done in two different ways: with respect 
to the energy shell center ('global average', localization parameter (39) (3d = f3 g ), 
and with respect to the centers of individual distributions ('local average', (3 d = (3j). 
Besides, the average < (3d > over (3 d values from the individual distributions was 
computed. 

In the ergodic case A = 3.6 both average distributions for EFs are fairly close to 
the SC law - a remarkable result, because that law was theoretically predicted for the 
other distribution, that of GSa. More precisely, the bulk ('cap') of the distributions 
are very close to the limiting SC (35), except in the vicinity of the SC singularities. 
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Numerical values of the localization parameter (f3 g = 1.08, (3 l = 0.94, < (3 d >= 0.99) 
are in a reasonable agreement with scaled (3d = 0.97 for A = 3.6, Eq.(39). 

As expected, the GS structure is similar: (3 g = 1.07, (3 l = 1.06, < (3d >= 0.98. 

For finite q all the distributions are bordered by the two symmetric steep tails 
which apparently fall down even faster than the simple exponential with a char- 
acteristic width ~ b. The physical mechanism of the tail formation is a specific 
quantum tunneling via intermediate BSs [30]. The asymptotic theory of the tails 
was developed in Ref. [30,42,46]. Surprisingly, it works reasonably well even near the 
SC borders. 

The structure of matrix c mn is completely different in the localized case A = 0.23. 
The EF local average shows clear evidence for exponential localization with (3 l = 0.24 
which is again close to scaled (3 d = 0.21 for A = 0.23. However, the global average 
reveals a nice SC (with tails) in spite of localization {(3 g = 0.98). It shows that, at 
average, the EFs homogeneously fill up the whole energy shell. In other words, their 
centers are randomly scattered over the shell. 

Unlike ergodic case the localized GS structure is quite different from that of 
EFs. Both averages now yield similar results which well fit the SC distribution 
(P g = 0.98, ft = 0.96 as compared with JJ g = 0.99, ft = 0.24 for EFs). So, GSa 
look extended, yet they are localized! This is immediately clear from the third 
average < (3 d > = 0.20. The explanation of this apparent paradox is that even 
though the GSa are extended over the shell they are sparse that is contain many 
'holes'. 

In analysis of the WBRM structure theoretical expression (38) for the ergodic 
localization length d e (the energy shell width) was used. In more realistic and 
complicated physical models this might impede the analysis. In this respect the new 
method for direct empirical evaluation of d e , and hence the important localization 
parameters (3d and A in Eq.(39), from both average distributions for GSa as well as 
from the global average for EFs [41] looks very promising. 

The physical interpretation of this structure based upon the underlying chaotic 
dynamics is the following. Spectral sparsity decreases the level density of the oper- 
ative EFs which is the main condition for quantum localization via decreasing the 
relaxation time scale (see, e.g., [7]). Yet, the initial diffusion and relaxation are still 
classical, similar to the ergodic case, which requires extended GSa. On the other 
hand, EFs are directly related to the steady-state density [11], both being solid 
because of a homogeneous diffusion during the statistical relaxation. 

This picture allows to conjecture that for a classically regular motion the EFs 
become also sparse, so that EF/GS symmetry is apparently restored. 
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5 Conclusion: pseudochaos and traditional 
statistical mechanics 



The quantum chaos is a particular but most important example of the new generic 
dynamical phenomenon - pseudochaos in almost periodic motion. The statistical 
properties of the discrete-spectrum motion is not a completely new subject of re- 
search, it goes back to the time of intensive studies in the mathematical foundations 
of statistical mechanics before the dynamical chaos was discovered or, better to say, 
was understood (see, e.g., Ref.[55]). This early stage of the theory as well as the 
whole TSM was equally applicable to both classical and quantum systems. For the 
problem of pseudochaos one of the most important rigorous results with far-reaching 
implications was the statistical independence of oscillations with incommensurate 
(linearly independent) frequencies cu n , such that the only solution of the resonance 
equation 

N 

™n ■ uj n = (41) 

n 

in integers is m n = for all n. This is a generic property of the real numbers. In 
other words, the resonant frequencies (41) form a set of zero Lebesgue measure. If we 
define now y n = cos (u n t) the statistical independence of y n means that trajectory 
y n {t) is ergodic in iV-cube \y n \ < 1. This is a consequence of ergodicity of the phase 
trajectory <p n (t) = u) n t mod 2n on a torus \(p n \ < ir. 

Statistical independence is the basic property of a set to which the probability 
theory is to be applied. Particularly, the sum of statistically independent quantities 

N 

X(t) = A n ■ COS (Uj n t + <f> n ) (42) 

n 

which is the motion with discrete spectrum, is a typical object of this theory. How- 
ever, the familiar statistical properties like Gaussian fluctuations, postulated (di- 
rectly or indirectly) in TSM, are reached in the thermodynamic limit N — > oo only 
[55]. In TSM this limit corresponds to infinite-dimensional models [6] which provide 
a very good approximation for macroscopic systems, both classical and quantal. 

What is really necessary for good statistical properties of almost periodic motion 
(42) is a big number of frequencies N u — > oo which makes the discrete spectrum 
continuous (in the limit). In TSM the latter condition is satisfied by setting jV w = 
iV — ■> oo. The same holds true for quantum fields which are infinite-dimensional. In 
the finite-dimensional quantum mechanics another mechanism, independent of N, 
works in the quasiclassical region Q 1. Indeed, if the quantum motion (42) (with 
if)(t) instead of x(t)) is determined by many (~ Q) eigenstates we can set N w = Q 
independent of N. The actual number of terms in expansion (42) depends, of course, 
on a particular state ip{t). For example, if it is just an eigenstate the sum reduces 
to a single term. This is reminiscent to some special peculiar trajectories of classical 
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chaotic motion whose total measure is zero. Similarly, in quantum mechanics N w ~ 
Q for most states if the system is classically chaotic. 

For a regular motion in the classical limit the quantity N u <C Q becomes consid- 
erably smaller. For example (see Ref. [14]), in the standard map = Q in ergodic 
case, JV W ~ k 2 in case of localization (both classically chaotic, K > 1) but only 
Nu, ~ k <C A; 2 <, Q for classically regular motion (if < 1). The quantum chaos-order 
transition is not as sharp as the classical one but the ratio N U (K > 1)/N W (K < 
1) ~ k — > oo increases with quantum parameter k. 

Thus, with respect to the mechanism of the quantum chaos we essentially come 
back from the ergodic theory to old TSM with exchange of the number of freedoms iV 
for quantum parameter Q. However, in quantum mechanics we are not interested, 
unlike TSM, in the limit Q — > oo which is simply the classical mechanics. Here, 
the central problem is the statistical properties for large but finite Q. This problem 
does not really exist in TSM describing macroscopic systems. In a finite-Q (or 
finite— N) pseudochaos we have to introduce the basic conception of time scale [11]. 
This allows for interpretation of quantum chaos as a new dynamical phenomenon, 
related but not identical at all to the classical dynamical chaos. Hence, the term 
pseudochaos emphasizing the difference from the asymptotic (in time) chaos in the 
ergodic theory. 

In my opinion, the fundamental importance of quantum chaos is precisely in 
that it reconciles the two apparently opposite regimes, regular and chaotic, in the 
general theory of dynamical systems. The studies in quantum chaos help to better 
understand the old mechanism for chaos in many-dimensional systems. Particularly, 
the existence of characteristic time scales similar to those in quantum systems was 
conjectured in Ref. [7]. 

Is pseudochaos a chaos? 

Until recently even the conception of classical dynamical chaos was rather incom- 
prehensible, especially for physicists. I know that some researchers actually observed 
dynamical chaos in numerical or laboratory experiments but... did their best to get 
rid of it as some artifact, noise or other interference! Now the situation in this field 
is upside down: many researchers insist that if an apparent chaos is not like that in 
the classical mechanics (and in the existing ergodic theory), then it is not a chaos at 
all. Hence, sharp disputes over the quantum chaos. The peculiarity of the current 
situation is that in most studies of the 'true' (classical) chaos the digital computer 
is used where only pseudochaos is possible that is one like in quantum (not classical) 
mechanics! 

Hopefully, this 'child disease' of quantum chaos will be over before long... 

Acknowledgment. The conception of quantum chaos presented above has been 
developed in a long-term collaboration with G. Casati, J. Ford, I. Guarneri, F.M. 
Izrailev and D.L. Shepelyansky. 
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